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O ; Abstract 
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A modified Gersch-Rodriguez formalism describing the leading Final State 
Effects in the high momentum transfer response of low concentration 3 He- 4 He 
mixtures is presented and discussed. The leading corrections to the Impulse 
Approximation are expressed in terms of the interatomic potentials and the 
semidiagonal two-body density matrices of both the mixture and its boson- 
boson approximation, in which the 3 He atoms are replaced by bosons of the 
same mass and at the same partial density. Numerical calculations of the Fi- 
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nal State Effects functions of He and He are finally presented and discussed 
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The dynamic structure function S(q, u>) is known to contain the maximum information 
a neutron scattering experiment can provide about the structure of a quantum fluid. Many 
theoreticalac! and experiment altTa work designed to study the Bose-Einstein condensate 
in superfluid 4 He has been done in the last thirty years, and much information regarding 
its quantum properties have been gathered so far. An alternative system in which Bose- 
Einstein condensation sets in is the isotopic 3 He- 4 He mixture, stable at low temperatures 
and concentrations, and presenting much richer static and dynamic properties due to the 
strong interaction between bosons ( 4 He) and fermions ( 3 He). However, 3 He atoms are high 
neutron absorbers and the experimental measurement of the high q dynamic structure func- 
tion of systems containing 3 He is technically difficult. For this reason, only one measurement 
of the high q response of the mixture has been reported up to the datetO, and the analysis of 
the data led to surprising conclusions such as a condensate fraction value no ~ 0.2, in clear 
disagreement with almost all theoretical predictions which are usually close to n « O.ltm. 
However, and despite the existent discrepancies between theory and experiment, no new 
measurements of the high momentum transfer response of the mixture have been performed 
yet, and hence the question of what this value really is still remains open. Consequently, the- 
oretical investigation of the dynamic structure function of the mixture and further analysis 
of the measured data is well justified. 

The theoretical analysis of the q — > oo response of the mixture is however slightly more 
involved that in the pure 4 He case due to the presence of a fermionic component. In fact, 
most of the formalisms that have been used in the past to describe Final State Effects 
(FSE) in pure 4 He rely, in one way or another, on the presence of long-range order in the 
the one-body density matrix pi(r). Long range order makes pi(r) reach a finite r — > oo value 
without changing sign, and thus allowing a cumulant expansion of the Fourier transform of 
S(q,u) around its high q limit pi(tq/m). This is not the case when dealing with fermionic 
components, and therefore none of these formalisms can be directly applied to analyze the 
high q response of 3 He without further modifications. In this work we address this problem 
and present a new formalism designed to describe FSE in the response of the 3 He- 4 He 
mixture, where a 3 He component is present but dilute. 

In the mixture and when the transferred momentum is high enough for the incoherent 
approximationOJa to hold, the dynamic structure function becomes 

S(q,u) = 04 (1 - x 3 ) S®(q,u) + * 3 x 3 S®{q,u) , (1) 

where x 3 is the 3 He concentration, and a 3 and 04 are the cross sections of the separate 
scattering processes (03 = 5.61 and a 4 = 1.34 in units of banish). Each of the individual 
responses are the Fourier transform of their corresponding density-density correlation factor 

1 N <* , 

S {a \q,t) = W J2 (e-^e^y^e-^*) (2) 

where a = 3 for 3 He and a = 4 for 4 He. In this expression, H stands form the Hamiltonian 
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corresponding to a mixture in which a 3 He and a 4 He components coexists and interact 
through pairwise local and central potentials V^ a '^(r). 

Actually, S^ (q, t) is the starting point from where a high q expansion of the response 
is performed, and this is done acting with the two translation operators exp(— iq ■ r^) and 
exp(zq ■ Tj) on the time evolution operator exp(-iHt) appearing in between. This leads to 

S («)( q , t ) = -U^' £ U^ + Lf) e -it H \ ^ (4) 

where 1/f = v^ a ' p.,, is the projection of the momentum of particle j along the direc- 
tion of the recoiling velocity v^ = q/m a and u)jf*> = q 2 /2m a may be understood as its 
kinetic energy when the momentum transferred to the system is high. Taking into ac- 
count the symmetries of the ground state wavefunction, introducing an identity in the form 
l = exp(— iLjt) exp(iLjt) and moving the time evolution operator to the left, the density- 
density correlation factor becomes 

S^\q,t) = e i ^ a)t (c {a \t)e itL ^\ , (5) 

where C^^s) is the operator carrying the effects of the FSE, and reads 

C^(t)=e- UH e u{H+L ^e- ltLl . (6) 



In the high q limit, v ( a > is large and t is short while their product remains of order 
unity. Therefore, one is allowed to write the response as a function of s = v^ a H instead of t, 
and expand it in inverse powers of the recoiling velocity v^ a \ Defining a new Hamiltonian 
H = Ejv {oL \ Eqs. (|) and @ become 

S<°0(g, s ) = eH 00 /* " (C^(s)e^ (a) ^) (7a) 

C w (s) = e -^ e «(*+*c«).pi) e -favM.pi ? (7b) 

where v^ is the unit vector pointing in the direction of \~( a \ Direct differentiation of (|7aD 
leads to the following first order differential equation satisfied by C^ a '(s) 

-^-C (a) (s) = C {a) (s) = % \H(s) - H] C {a) (s) (8) 
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which must be solved together with the initial condition C^(0) = 1, the latter being implicit 
from the definition of C^ a '(s) in Eq. (|7b|) . In Eq. (§), W(s) is an effective Hamiltonian that 
results from the action on TC of the different exponential operators appearing in C^ a '(s) 

H(S) = e -isH e i S (H+^. pi ) He -i S (H+v("1. pi ) e is<H (9) 

Different representations of the density-density correlation factor can be obtained from 
the different forms in which the solution of Eq. (^) can be organized. One form that has 
proved to be particularly useful in the past is the cumulant onell3o, which captives the most 
salient features of the response and simplifies in the q — > oo limit 



C {a \s) = e ^A ( S ) e 6 2 a A 1 ( S ) e 4A 2 ( S ) . . . = JJ exp [ e "A n ( 5 )] 



(10) 
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Assuming this analytical form, the first s-derivative of C^ a \s) can be easily calculated and 
plugged into Eq. (|[). Then equating each order in e™ in both sides of the equation one 
finds a hierarchy of coupled differential equations that can be solved for every A n once the 
previous AA^ . . . , A n _i are known. The first terms in the chain can be readily obtained 
and lead to e„ 



C (a \s) = exp [ze a J* (H (y) - H) dy 

H 2 (s')ds' - 
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Alternatively, C(s) can be casted as an additive series setting 

e^ A ^ = l + e n a T n (s) 



ds'ds' 1 



;n) 
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where T n (s) is to the lowest order independent of e a . In this case one finally arrives at the 
general solution 



C^(s) = [l + e a r 1 (s) 



i + eir 2 (s) 



i + eir 3 (s) 



Je a f°(H (s>)-H)ds> , ( JelJ'H^s'W 
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Equations ( |TTT ) and (|13D are different expressions of the same solution. Both of them can 
be used in Eq. ([Ta|) to generate a representation of the density-density correlation factor in 
which each order in e a is separately displayed. As a matter of fact, both forms coincide in 
the e a — > limit, where only the zero order term contributes. In this case C^ a '(s) = 1 and 
S^ (q, s) reduces to the Impulse Approximation 
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which is known to dominate the q — > oo behavior of the response. At high but finite q, e a 
is small and thus the leading corrections to the IA are already captived by the first term in 
Eq. (|T3"D, and so all other terms may be discarded. This leads to the following approximation 
for the response 



S[ a \q jS )=e* 



is4 a) /v^ I ie a f°(H (s')-H)ds> ts-pi 



(15) 



The difference Hq(s') — H in the argument of the exponential can be written in terms of 
the interatomic potentials alone. In the general case considered here in which two different 
species coexist and interact through pairwise local and central forces, the argument of the 
exponential reduces to 



H a (s) - H 



D «s-pi 



He 



ispi 
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E E^K' + s)-^^^,))^ E E^ (a/3 Ws) 



(16) 
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where in the last line use has been made of the definition of the potential difference 

V^(r„ + rO-V (Qffl (r«). (17) 
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Moreover, the leftmost exponential operator inside the expectation value in Eq. (|T5|) is 
diagonal in configuration space, and therefore S\ (q, s) reduces to an integral over the 
semidiagonal iV-body density matrix and the interatomic potentials 



S[ a \q,s) 



N 3 \N A \ 



e *4"V«M /dr^ JV (r 1 ,r 2 ,...,r JV ;r 1 H-s) 



x exp 



v (a) j 



-EE/ ds> AvW\r Xj ,s>) 
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At zero temperature p^ is defined in terms of the ground state wavefunction ty 

p N (r u r 2 , . . . , r N ; r[) = N 3 \N 4 \%(r u r 2 , . . . , r N )^f(r[ 7 r 2 , . . . , r N ) , (19) 

where the N 3 coordinates corresponding to 3 He particles represent both position and spin 
coordinates. 

Eq. ( |18D is simple compared to the expression of the exact response, but it is still very 
difficult to evaluate due to the iV-body quantities entering on it. However, it is in the 
appropriate form for a Gersch-Rodriguez cumulant expansion!!^. While the original Gersch- 
Rodriguez formalism deals with expectation values of time-ordered integrals of operators, 
the approximations made so far have brought the response to the Static Background Ap- 
proximation. In this limit, the scattering time is assumed to be so short and the momentum 
transfer so high that in essence only the particle being struck by the incoming neutron moves 
in the process, and thus all other particles in the background are assumed to be frozen at 
their positions. Under these circumstances, the series simplifies and reduces to the general 
expansion rule 
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(20) 



where 0(i"i, r^), /(ri, . . . , r^; r' : ) and 9 m (s) are arbitrary functions, the latter also possibly 
depending on particle coordinates r 1; r' 1; r 2 , . . . , r N . W and u n are the coefficients of the 
expansion and up to n = 1 read 



Wa 
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(r 1 ,r; ) + /^/(r 1 ,r 2 ,...,r„;r; 
-— E dr N f(r u ..., r N ;r[ 

VVn m J 



exp 



ds'9 m (s' 



(21a) 
(21b) 



while expression for higher order terms can be derived following the steps described in the 
literature^!. Notice that due to the logarithmic nature of the expansion, this relation holds 
only if Wq has no zeros. In the current case, Wq is related to the one-body density matrix 
of the component whose density-density correlation factor is being analyzed, and whereas 
p\ (s) is everywhere positive and fulfills the required condition, p[ (s) presents a complex 
nodal structure that makes the applicability of the expansion rely on the appropriate choice 
of 0(ri,r'i). 

When expressions (|20|), ( |21a|) and (|21b| ) are used to compute the density-density corre- 
lation factor of the 4 He component of the mixture, one arrives at the natural extension of 
the original Gersch-Rodriguez result. Setting (frfji,^) = 0, one gets 
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W = / dr N p N (r 1 ,r 2 ,...,r N ;r l + s) = —p ( i\s) 

1\3\1\4\ J pi 



(22) 



while uj\ reads 

1 p 4 
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N 3 \N 4 \ p W {s ) 



£ E [dr N p N (r 1 ,...,r N ;r 1 + s) 1 - exp ( * f ds' AV^\v Xj , s> 



5=3,4 j'e/3 



(4) 
Pi 



drp^ ' (r, 0; r + s) 



1 «!>(^/ o ds' AV^(r,s>) 
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Pi 



drp 



(4,3), 



r, 0;r + s) 



eXP (^)/o S ^ AV(43)(r ' S ' ; 



(23) 



in terms of the (4, 4) and (4, 3) components of the semi-diagonal two-body density matrix 
of the mixture 



(«,/3)/ M 

p\ ; (ri,r 2 ;r ly 



N a (Np - 5^] 



Jdr 3 --- dr N %(r 1: r 2 , . . . , rjv)^o r i> r 2, • • • , t n ) 



Jdv N ¥<n, r 2 , . . . , r N ) 



(24) 



At this level, therefore, S^(g, s) is predicted to be the algebraic product of the IA and 
the FSE broadening function 
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S[ A \ q ,s) = S\%,s)R (4 \ q ,s) 



(25) 



where 



Pi 



(26) 



and 



1 



L P?(s) 



J drpf A) (r, 0; r + s) 1 - exp l-^ f ds'A V (44) (r, s'; 



Pi ( S , 



f rfrp^ 4,3) (r,0;r + s) 



l-exp(— - / ds' AV^\r,s') 
v W JO 



(27) 



These results are formally equal to the Gersch-Rodriguez ones used to compute the high q 
response of pure 4 He, the only difference being the presence of p 2 ' which is the contribution 
coming from the interaction of the 4 He atoms with the 3 He atoms in the mixture. As a 
matter of fact, in the zero 3 He concentration limit this last term cancels and expression (J27D 
coincides exactly with the one reported in ref.^, thus stressing that the former is the natural 
extension of the latter to the mixture where two different isotopes coexist. 

Unfortunately, no cumulant expansion of the 3 He density-density correlation factor can 
be performed as described above. This is because the zero order cumulant, which is pro- 
portional to the 3 He one-body density matrix, has infinitely many nodes. This problem can 
be bypassed recalling that in the high momentum transfer limit the most relevant processes 
in the scattering are those taking place at short distances, where dynamical correlations 
dominate over statistical ones. This means that in the q — > oo limit, and disregarding the 
effects on the Bose-Einstein condensate, FSE in boson and fermion systems should look 
like similar, and that therefore the contribution of the fermion statistics to the FSE can be 
introduced as a small correction to the effect produced by the dynamical correlations, which 
are entirely taken into account by a purely bosonic FSE function. But on the other hand, 
the IA is known to substantially depend on the statistics obeyed by the system. The s rep- 
resentation of the IA is proportional to the one-body density matrix, so it seems that this is 
the only quantity that really requires a proper treatment of the statistics. One can consider, 
therefore, a factorization of the iV-body density matrix of the mixture in the following form 



pjv(ri, r 2 , . . . , r N ; ri) = p\ >(r iv ) 



(3), 



1 



pfOn 



-pN(r 1 ,r 2 ,...,r N ;r[) 



+ 



t /x Pi (rii') b, m 

Pi(rii') 



(28) 



where pf and p^ are the one- and the iV-body density matrices of a fictitious mixture 
(henceforth referred to as the boson-boson mixture) in which the 3 He atoms are replaced by 
bosons of the same mass and at the same partial density. Notice that this replacement does 
not introduce any singularity in p^ because pf (r) corresponds to a bosonic phase and thus 
has no nodes. 

With this prescription, the 3 He response of the mixture becomes the sum of two terms 



y(3) 



S[ 3 \q,s) = S%\q,s) + ASW(q,s), 



(29) 



where 



S { B } (q,s) 



e ^>/„(3)l p (3) 
P3 



(30) 
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P3 



N 3 \N A 
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N 
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Pf« 



p N {r l ,r 2 ,...,r N ;r 1 + s) exp 
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and 



AS (3) {q,s) = e 



isuf ] /v^\ 
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N*\NA 
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bi \ Pn( t i, r 2, • • • , r N ; n + s 



p N (r 1 ,r 2 ,...,r N ;r 1 + s) (31) 

£ WVa^^s') 

> n a a i- a " U 



exp 






0=3,4 je/3' 



In this approximation, Eq. fl3~0|) describes that part of the response that can be written 
as the product of the IA, which enters through p\ \s), with a bosonic FSE function, that 
is given by the term in square brackets. It is important to notice, however, that the IA 
is exactly accounted in this term because the exact 3 He one-body density matrix has been 
previously factorized; while the FSE broadening function is evaluated in the boson-boson 
approximation. In contrast, AS^(q, s) carries the contribution of all the statistical correla- 
tions between particles 2, 3, . . . , A3 of 3 He that have not been taken into account in p\ (s). 
This last terms is expected to introduce small corrections to S^ (q, s) that only appear when 
the response is computed beyond the IA level. 

As before, Eqs. ( J30|) and (J3ID are difficult to handle due to the presence of p^v and p^. 
Nevertheless, the FSE function in ( p0|) can be worked out just as in the 4 He case due to the 
bosonic nature of the functions entering on it. Hence one finds 



7(3) 



S^\q,s) = S ( I 3 >(q,s)R^(q,s) 



where 



S < I %,s)=e* 



4 3) M 3 >± p (3) (s) 

P3 



(32) 



(33) 



and 

i? (3) (g,s) = exp 



P?{s) 



drp^ (r, 0; r + s) 



1 — exp 



^0) 



P?(s) 



drp^ (r, 0; r + s) 



ds' 'AV (33) (r,s' 

expf-^r f 8 ds' AF( 34 )(r,s' 
\v( 6 > Jo 



(34) 



The 3 He additive term AS^ 3 '(q,s) can not be handled in the same way because in this 
case the zero order cumulant vanishes 



W 



1 



A 3 !A 4 ! 



dr 



N 



p N (r 1 ,r 2 ,...,r 1 + s) 



P?(s)n 



Pf(s) 



p 7V (r 1 ,r 2 ,...,r 1 + s) 



(35) 



thus violating the condition imposed on W . 

The problem of finding a useful prescription for AS^(q, s) can be solved inspecting the 
structure of the two-body density matrices that would enter in the lowest order terms of 
a cumulant expansion, as they carry the leading contributions to the FSE in the q — > 00 
limit. When Eq. fl28D is integrated over all particle coordinates but i - ! and r 2 , an equivalent 
factorization of p 2 is found. The difference between the 3 He semidiagonal two-body 
density matrices of the mixture and their boson-boson approximation 
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(3,a) / /\ 

P2 (ri,r 2 ;r ly 



.(3) 



Pi V ) „(3,a) 

PiViv 



P2 



( r i 5 r 2;ri, 



(36) 



can be analyzed in the framework of the HNC/FHNC equations starting from a variational 
model of the ground state wave function. Careful inspection of the diagrams entering in fl36|) 
reveals that Ap 2 (ri, r 2 ; r'i) may be factorized as followstZrc2 

Papf(r n ,)G^\v Xl T 2] v[) - p a p 1D {r xv )F^ a \ 



Apf a \ 



ri,r 2 ;r 1 ; 



r i; r 2; i"i, 



(37) 



where Pld(?~ii') is an auxiliary function that adds the contribution of all those diagrams 
linking points 1 and V that are not connected to point 2 and with no statistical lines starting 
or ending at points 1 and 1'. Function G ( - 3 ' a \r 1 , r' x ; r 2 ) and F^ 3 '°^(ri, r[; r 2 ) are the sum of all 
other diagrams not contributing to p[ (s) that contain dynamical and statistical lines linking 
the external points. Indeed, it can be seen from its diagrammatic structure that Pid(^ii') 
shares many common features with the one-body density matrix of a purely bosonic liquid, 
as for instance it is always positive and its large r value approaches a constant that would 
be ascribed to some sort of condensate fraction value, although it has no physical meaning 
in this case. As a matter of fact, numerical calculations show that at low 3 He concentrations 
Pid{t\v) and pf (riy) are nearly identical, and so that both P\d{t\v) and pf (riy) satisfy the 
expansion condition and that either of them can be used as the basic function 0(ri, r' x ) upon 
which the cumulant expansion is being built. Choosing Pid(^h') as the starting function, 
AS^ 3 \q, s) in Eq. (|3~TD can be brought to a form suitable for expansion purposes by adding 
and subtracting the former to the latter 



AS®(q,s) = e i ™ W ' v&) 



1 

P3 



Pid(s) + N * N j J dr N ( pjffa, r 2 , . . . , r^; r x + s) 



P ( l\s) (3,3)B, 

Pid{s) 



r 2 ,...,r^r 1 + s))exp(^- £ £ f ds' AV^^s' 



/3=3,4 jeP 



-e^ 3) ^ (3) Ip 1D ( s) 

P3 



(38) 



as now the term inside the square brackets admits a cumulant expansion. Up to the first 
order this leads to 



AS^(q,s) = e^ ) ^ ) ^—pf\s) 



exp 



Pid{s) 
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(39) 
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f drAp 2 3 ' 3) (r, 0; r + s) 1 - exp ( -^ f ds'AV {33) (r, s 



Pid(A 



f drAp { 2 3 ' 3) (r,0;r + s) 



1 (>x l>(^y/ ds' AV (33 \r,s] 
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The results in Eqs. @, (|33), (§§, (|3J) and ((5|) constitute the prediction of the 3 He re- 
sponse of the mixture in the present formalism. As before, the FSE functions are built upon 
the interatomic potentials and the semidiagonal two-body density matrices, even though 



the latter must now be computed for both the real mixture and its boson-boson approxi- 
mation. As in the 4 He case, these results resemble the original ones derived by Gersch and 
Rodriguez and form in fact their extension to the mixture when the problem of finding a 
useful prescription for the 3 He response is addressed. 

The 4 He and 3 He FSE functions can be evaluated once one has a suitable description 
of the required ground state two-body density matrices. These can be computed starting 
from a variational Jastrow wavefunction in the framework of the HNC/FHNC equations 
for the mixture, generalizing the formalism developed in Refs.E-j'E3. In the simplest scheme 
one disregards the contribution of the Abe diagrams and works in the Average Correlation 
Approximation (ACA), in which the (4,4), (4,3) and (3,3) correlation factors are assumed 
to be equal. In this limit, differences between the isotopes are left to the effects derived from 
their different mass and statistics. Despite the simplifications, the wavefunction generated 
in this way still captives the essential features of the mixture at T = 0. 

In the particular case of the 3 He- 4 He mixture, the interaction does not distiguish between 
isotopes and so all three potentials JA a, ^(r) have been taken to be equal to the HFDHE2 
Aziz potential^. The computed FSE functions of the 23 = 0.095 3 He concentration mixture 
at equilibrium density p = 0.3554 cr -3 are shown in figure (1). The 4 He ( 3 He) FSE function 
has been Fourier transformed from s to the I4 (Y3) West scaling, where Y a = m a uj/q — q/2. 
The solid and dashed lines on the left show R^ (q, Y4) compared to the Gersch-Rodriguez 
FSE function of pure 4 HetM at saturation density (p = 0.365a™ 3 ). The solid, dashed 
and dotted lines on the right depict R^(q,Y 3 ), AS^(q, Y3) and the 3 He Compton profile 
J^iYs) — (l/'m'3)S^(q,Y 3 ) in the IA. As it can be seen and despite the different densities, 
R^(q, Y4) is quite similar to the pure 4 He FSE function, the main differences being present 
at the tails. On the other hand, R^ (q, F 3 ) is close to both functions, even though the peak is 
slightly higher. This effect is mostly due to the low partial density of the 3 He component at 
the mixture p^ = x^p = 0.0338 a™ 3 , compared with the partial density of the 4 He component 
P4 = (1 — x^)p = 0.3216 cr -3 and the equilibrium density of pure 4 He. Moreover, the 3 He 
additive term is rather small compared with the IA prediction, thus indicating that at high 
q statistical corrections to the purely bosonic FSE in the 3 He peak are present but play a 
much less significant role compared to the effect produced by R^(q,Y 3 ). 

In summary, we have presented a new formulation of FSE for 3 He- 4 He mixtures, where 
the presence of a fermionic component forbids a straightforward application of most of 
the existing FSE theories used in the analysis of the high momentum transfer response of 
pure 4 He. The formalism is inspired in the theoretical analysis carried out by Gersch and 
Rodriguez, and actually reduces to it in the zero 3 He concentration limit. The resulting 
expressions are expected to accurately describe the effect of FSE in the response of the 
mixture. Furthermore, numerical calculations reveal that the 4 He FSE function is similar to 
the pure 4 He one. In this way, FSE corrections to the 4 He peak in the mixture and in the 
pure phase are expected to be similar. Statistical effects do not appreciably modify a picture 
in which the FSE of the 3 He component are computed in the boson-boson approximation. 
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FIGURES 

FIG. 1. Left: R( 4 '(q, Y4) compared to the pure 4 He FSE function of Gersch and Rodriguiez 
(solid and dashed lines). Right: R^ 3 '(q, Y3) and AS^(q,Y^) compared to the IA prediction for the 
3 He peak in the mixture (solid, dashed and dotted lines, respectively). 
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